Interaction-induced anomalous transport behavior in one dimensional optical lattice 



Zi Cai 1 , Lei Wang 1 ^ X. C. Xie 2 ' 1 , and Yupeng Wang 1 
1 Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, 
Chinese Academy of Sciences, Beijing 100080, P. R. China and 
2 Department of Physics, Oklahoma State University, Stillwater, Oklahoma 74078, USA 

(Dated: November 18, 2009) 

The non-equilibrium dynamics of spin impurity atoms in a strongly interacting one-dimensional 
(ID) Bose gas under the gravity field is studied. We show that due to the non-equilibrium prepa- 
ration of the initial state as well as the interaction between the impurity atoms and other bosons, a 
counterintuitive phenomenon may emerge: the impurity atoms could propagate upwards automat- 
ically in the gravity field . The effects of the strength of interaction, the gradient of the gravity 
field, as well as the different configurations of the initial state are investigated by studying the 
time-dependent evolution of the ID strongly interacting bosonic system using time-evolving block 
decimation (TEBD) method. A profound connection between this counterintuitive phenomenon and 
the repulsive bound pair is also revealed. 
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Recently, ultracold atoms in optical lattice have pro- 
vided a perfect platform for simulating quantum many- 
body models in condensed matter physics [TJ [2]. What's 
more, the uniqueness of cold atomic system, such as the 
low dissipation rate as well as long coherence times, has 
opened exciting possibilities for studying non-equilibrium 
quantum dynamics of many-body systems. It is known 
that, due to the energy dissipation between the system 
and the environment, the properties of the many-body 
system in solid state physics is mainly determined by its 
ground state as well as its low-energy excited state. In 
the ultracold atomic systems, however, the extreme low 
dissipation rates there guarantees the conservation of the 
system energy and total particle number in a relatively 
long time. Therefore, not only the ground state but also 
the high energy excited state may contribute to the non- 
equilibrium dynamics of the many-body systems, which 
may exhibit novel phenomena [3j EJ 13 E 13 IE] • 

Repulsively bound atom pairs is one of the most in- 
teresting and novel phenomena emerging from the non- 
equilibrium many-body physics in optical lattice [9j [TO] . 
It is shown that though repulsive force separates par- 
ticles in free space, under a periodic potential and in 
the absence of dissipation, a bound atom pairs could 
be stabilized by the strongly repulsive interaction due 
to the conservation of energy. This unconventional phe- 
nomenon provides a typical example of non-equilibrium 
physics determined by the high-energy excited state of 
the Hamiltonian rather than the ground state, and has no 
analogue in traditional condensed matter systems due to 
the rapid dissipation. In this Letter, we provide another 
example of these exotic non-equilibrium phenomena by 
studying the propagation of spin impurity atoms through 
a strongly interacting one-dimensional (ID) Bose gas un- 
der an external potential which decreases linearly along 
a particular direction. The ground state and the quench 
dynamics of the Bose-Hubbard model under linear exter- 
nal potential have been discussed previously pT| \12\ . We 



show that if the initial state is prepared far away from 
the equilibrium ground state, the spin impurity atoms 
could propagate towards the direction the external lin- 
ear potential increases, which means a particle with ini- 
tial momentum zero moves upwards in a gravity field in 
our case. This unconventional phenomenon, similar to 
the repulsive bound atom pairs, is a result of the conser- 
vation of the energy. We would also discuss the profound 
connection between these two non-equilibrium phenom- 
ena. 

First, we propose the experimental set-up of our 
non-equilibrium system. It is motivated by a recent 
experiment [13] , where the quantum transport of spin 
impurity atoms through a Tonks-Giradeau gas is stud- 
ied. The departure point of our experimental set-up is 
loading the 87 Rb atoms with the hyperfine ground state 
\F = l,77ip = —1) into a ID optical lattice, which is 
along z-direction. We assume that the interaction be- 
tween these 87 Rb atoms is very large and the filling factor 
is one, thus the ground state of the system is the Mott 
state [14]. An additional harmonic magnetic trap is im- 
plemented to provide a vertical confinement for the atoms 
and prevent the atoms escaping from the optical lattice 
due to the gravity. The harmonic magnetic trap cou- 
ples with the atoms via the magnetic dipole interaction, 
therefore the vertical confinement is purely magnetic. 

To prepare the initial state, one pumps the Mott 
ground state into a non-equilibrium state by applying a 
pulse of radio frequency resonant with the energy gap 
between the two hyperfine spin states of 87 Rb atom: 
\F = l,rriF = —1) and \F = l^mp = 0) and inducing 
a transition \F = l^mp = —1) —> \F = l^rrtp = 0)[T5]. 
Below we denote the \F = l,mj7 = —1) as trap atoms 
and \F — l,mf = 0) as impurity atom. The spa- 
tial width and position of the pulse could be controlled 
experimentally [13]. In our case, we constraint the pulse 
to only induce the transition on one site and produces one 
impurity, as shown in Fig.l (general case containing more 
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FIG. 1: (Color online) (a). Experimental setup. The ini- 
tial state is prepared by applying rf pulse on single site, 
making boson atom transitions from \F — 1,?tif — 1) to 
\F = l,rriF = 0) level. Bosonic atom in the latter level expe- 
riences gravity potential, (b).-(d). Evolution of density dis- 
tribution of the spin impurity, for G/t = 0.5, U/t = 0,2,— 2 
respectively. 
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FIG. 2: (Color online) The COM of the spin impurity for (a). 
G/t = 0.5 with different interaction strengths, the U/t values 
are indicated at right axis. (b). U/t = 1 with different gravity 
fields, which could be attained by changing the lattice spacing 
or choosing different type of trapped atoms. 



than one impurities will be discussed below). Notice the 
magnetic moment of the impurity atom \F = 1, nip = 0) 
is zero thus it is not confined by the harmonic mag- 
netic trap and experiences an external linear potential 
produced by gravity. The experimental set-up proposed 
here is different from that in Ref.[l3] in the way that in 
our case, the bosons are loaded in a ID optical lattice 
along the z-direction, rather than the Tonks-Girardeau 
gas loading in a continuous ID space. Without the con- 
straint of the periodic potential along the z-direction, an 
impurity atom in a gravity field would be accelerated 
downwards and finally escape from the system, as shown 
in Ref . [T3] . In our case, however, the dynamics of the 
impurity atom becomes complex and nontrivial due to 
the optical lattice structure as well as the interaction be- 
tween the impurity atom and the trap atoms. Without 
the interaction, the single particle dynamics is known as 
Bloch oscillation [16], which has been observed in ultra- 
cold atomic systems [T71 H8]. 

The Hamiltonian in this system can be described by a 
ID two-component Bose-Hubbard model in a linear ex- 
ternal potential: 

i i i 

(i) 

where a^(^) is the annihilate operator for the impurity 
(trap) atoms. U a b — U is the interspecies on-site in- 
teraction strength. We assume the intraspecies on-site 
interaction U aa and Ubb are repulsive much larger than 
U a b, thus the boson could be considered as hard-core for 
the bosons with same species. G = mgao is the poten- 



tial gradient of the gravity, where m is the mass of 87 Rb 
atom and ao is the lattice constant of the optical lattice. 
The energy scale of G could be estimate as jiK, which 
is in the same order of magnitude of U and t. We as- 
sume that the energy scale of the gap between the s-band 
and p-band is much larger than that of the parameters 
(t,U,G) in our Hamiltonian. (1), therefore all the atoms 
are confined within the s-band and no interband physics 
involved. As we analyzed above, only the impurity atoms 
experience the gravity potential while the trap atoms are 
still confined by the magnetic trap. The non-equilibrium 
initial state can be represented as: |Tt , ''!' , 'tt)? 
where T (!) denotes the trap (impurity) atoms. The im- 
purity atom is initially localized in the center of the op- 
tical lattice, while all the trap atoms are still in the Mott 
phase. Obviously, this initial state is not a eigenstate of 
the Hamiltonian. (1) thus it would evolve. The time evo- 
lution of the system is controlled by the parameters in 
Hamiltonian. (1). Below we would study the dynamics of 
the impurity atom and show how it depends on the pa- 
rameters (U/t, G/t) in Hamiltonian. (1) using the TEBD 
method [19]. Total number conservation is used to reduce 
the computational effort. A straightforward treatment of 
the two component HCB needs local dimension d = 4, 
to reduce the computational effort, unfold technique is 
used and thus the local dimension is reduced to d = 2. 
The pay off is that the chain length increases to 2L and 
next-nearest neighbor interactions are introduced. Swap 
technique [I9j [20] is used to deal with next-nearest neigh- 
bor interaction within TEBD method. Both the unfolded 
(d = 2) and folded (d = 4) algorithm are implemented 
and to cross check the results. In this work we deal with 
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chain length L = 33, and the spin impurity was create 
at the 17th site. Finite size effect has little effect to the 
conclusion since the Bloch oscillation physics confines the 
atom near its original position. In the course of real time 
evolution we take the truncation dimension x = 80 and 
time step At = 0.05, the convergence is checked by tak- 
ing larger x- 

The time evolution of density distribution of the impu- 
rity atom is shown in Fig.l,(b)-(d). The dependence of 
the motion of impurity atom's center of mass (COM) on 
the interaction (U) and potential gradient (G) is shown in 
Fig. (2). If there is no interaction (U=0), the COM of the 
impurity atom does not move at all. This phenomenon 
could be understood from the dynamics of a single par- 
ticle under a combination of periodic and a linear ex- 
ternal potential. The presence of the external potential 
breaks the translational symmetry, render momentum k 
not a good quantum number. However, by choosing the 
temporal gauge (the scalar potential = 0), the exter- 
nal potential is taken into account by a time-dependent 
gauge field A(t) = Gt and the momentum k is replaced by 
k — A(t). The noninter acting Hamiltonian can be rewrit- 
ten in the momentum space as : H(t) = ^ k e k _A{t) a \ a k, 
where k — A(t) is good quantum number. Suppose ini- 
tially the particle locates on the point x = 0, the ini- 
tial state can be represented in the momentum space as: 
\x — 0)t=o = j- Y^k W) (e tkx =l because x = 0), where \k) 
is the single particle eigenstate in the momentum space 
and the sum is over the first Brillouin zone. In the pro- 
cess of time evolution, all the \k) modes experience Bloch 
oscillation from different initial momentum. The local- 
ization of the COM is actually a consequence of the su- 
perposition of all the momentum modes \k) in our initial 
state, the Bloch oscillation in \k) and | — k) modes cancel 
with each other and the net effect is zero. 

The single particle picture does not work when U ^ 0, 
thus we deal with the time evolution by TEBD. If U > 0, 
the COM of the impurity atom would move downwards 
until it reaches a steady quasi-equilibrium position, while 
interesting phenomenon emerges for U < 0. The COM 
of the impurity atom would move upwards though it is 
experiencing a gravity potential. We also find that the 
dependence of the motion of COM on U is not monotone. 
Notice in both case of large U and small U, the COM 
is always constrained around the original position. The 
maximum deviation of the COM away from original point 
occurs when U/t = ±2. The dependence of the dynamic 
of the COM on the gravity gradient G is also nontrivial. 
We notice that for larger G, the COM of the impurity 
would drop to a higher a final quasi-equilibrium position. 
The period of the oscillation is approximately the Bloch 
period 2tt/G. 

To understand the origin of the counter-intuitive phe- 
nomenon, we perform a Jordan- Wigner transformation to 
map our system into a interacting fermionic system (con- 
sidering the hard-core nature for the intraspecies boson) . 
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FIG. 3: (Color online) (a). Variance of the COM of spin 
down boson, (b). Total doublon numbers in the lattice for 
G/t — 0.5 and different interacting strength, U/t values are 
indicated at right axis. 

Then we use the particle-hole transformation on the trap 
atoms, thus the problem could be translated into a two- 
particle problem due to the half-filling condition in our 
initial state. The initial configuration in Fig[TJa) could 
be considered as two particles (an impurity atom and a 
trap hole) located on the same site (doublon) of the lat- 
tice. Below we mainly focus on the attractive condition 
(U < 0) in original model, which could be mapped into a 
two-particle problem with replusive interaction (U > 0) 
under the particle-hole transformation. If U ^> t, the 
dynamics of this two-particle problem reminds us of the 
physics in the the repulsively bound atom pairs [9j fTU] . 
The initial state is a doublon state with a high energy U. 
If the two particles are separated, the loss of the interac- 
tion energy (U) should be compensated by the increasing 
of the kinetic energy of the single particle to preserve the 
total energy of the system. However, in the periodic lat- 
tice structure, the kinetic energy of a single particle has 
a maximal value (2t). Therefore, when U ^> t, a doublon 
would be stabilized by the strong repulsive interaction, 
flgQc). In our case, interesting thing happens for smaller 
U, where the gravity field play a key role. It is possible 
that the loss of the interaction energy could be compen- 
sated by the increasing of the gravity potential energy of 
the impurity atom to preserve the total energy, as shown 
in figQb). The energy transfer between the interaction 
energy and gravity potential energy results in the anoma- 
lous dynamics of the impurity atom. 

To verify this qualitative picture, we again use TEBD 
to study the two particle dynamics from the given initial 
doublon state. The result is shown in Fig(3j We introduce 
the variance of the COM of the impurity atom, defined 
as Var = i 2 n^} — ((X^ n i) 2 )i to characterize the 
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FIG. 4: (Color online) (a). Sketch phase diagram of the final 
quasi-equilibrium state of the spin impurity: The impurity 
atom moving (I) upwards or (II) downwards for intermedi- 
ate attractive or repulsive interaction. Trapped exciton state 
induced by (III) strong attractive (IV) repulsive interaction, 
(b). Schematic view of impurity atom moving upwards to 
preserve the total energy, (c). Exciton be trapped by strong 
attractive interaction where \U\ ^> G. 



width of the wave packet, and the total doubly occupancy 
number D = n^n\) to study the fate of the initial 
doublon state, where n\ denotes the number of the trap 
hole on site i. The variance decrease as the increase of the 
interacting strength, although it shows oscillation which 
is given by the Bloch oscillation. This means that larger 
interacting strength squeeze the wave packet harder since 
dispersion of it lose more energy. For very large U (com- 
pared to G and t), Var is approaching to and D to 1, 
which means the doublon state is stabilized by the strong 
repulsive interaction. In addition, the doublon is local- 
ized on its original position since motion of it towards any 
direction violate the energy conservation in the gravity 
field. For smaller C7, interacting energy could be released 
to the energy continuum by pumping the impurity atom 
upwards to gain potential energy. This picture is further 
confirmed by the rapid decreasing of the total double 
occupancy number, which means the doublon dissolve 
shortly for small interacting strength. The COM and 
its variance of the impurity atom could be measured by 
time resolved tomographic measurement of the impurity 
density distribution [T3] . 

Translate back into spin impurity language, the im- 
purity atom and a trap hole bind together to form an 
exciton (which is the doublon in the two-particle pic- 
ture) for both large attractive and repulsive interaction. 
For the noninteracting case, the exciton dissolves imme- 
diately and Bloch oscillation physics dominates. For the 
intermediate interacting strength, the impurity atom will 
move upwards or downwards depending on the sign of the 
interaction. The sketch phase diagram is shown in Figj4] 



Notice that there is crossover instead of phase transition 
between different regions. The general principle guiding 
the time evolution of the strongly interacting systems 
is that the system always evolves toward a maximum 
entropy state with the constraint of the conservation of 
good quantum numbers (such as energy and total par- 
ticle number) [2TJ [22] . So the competition between the 
energy and entropy play a key role in determining the 
properties in the final quasi-equilibrium steady state. 

Finally, we would discuss the dynamic of the impu- 
rity atom evolving from two different initial configura- 
tions, we focus on the attractive condition to show the 
dependence of the anomalous phenomenon on the initial 
configuration we choose. The first case is that the fill- 
ing factor in the initial state is not one, which means 
the impurity would propagate in an environment of the 
quasi-superfluid of the trap atoms. The attractive inter- 
action between the impurity and the BEC atoms leads 
to a polaron[23], and its transport properties have been 
studied previously [24 . We choose the initial state as a 
quasi-superfluid with 31 hard-core bosons filled in the Id 
lattice with 33 sites. We can find that the COM of the 
impurity atom does not move upwards, different from 
that in the Mott initial state. This phenomenon can also 
be understood by the particle-hole transformation on the 
trap atoms. In this case, the trap hole is not localized in 
the initial configuration due to the quasi-superfluidity of 
the trap atoms, which means there is no high-energy dou- 
blon initial state. The energy transfer picture is no longer 
available, thus the anomalous transport phenomenon dis- 
appears. Actually, we can always prepare a Mott plateau 
region in lattice by adding an external harmonic trap po- 
tential, even the filling factor is not one, and the "Nega- 
tive Mass" phenomena will still appear. The other con- 
figuration is to apply a pulse of radio frequency on more 
than one site and produce many impurity atoms in the 
initial configuration. The initial state is chosen as five 
impurity atoms (from site 15 to site 19) located in the 
middle of the lattice while all other site is filled with 
trap atoms. We can find in this case the COM can still 
move upwards in the attractive interaction and the main 
result of this paper is not qualitatively changed by this 
multi-impurity effect. 
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